Thermal Analysis of the Pioneer Anomaly: 
A method to estimate radiative momentum transfer 



OO 
O 

o 

(N 
•*-> 

o 
O 
o 

(N 



X3 

Of 
i 

<L> ■ 
O ■ 

Oh! 

C/3 . 

o ■ 

>V 

^ : 

Ph. 



(N 

> 

o 
p 

o 

OO 

o 



X 



0. Bertolami 1 ' 3 , F. Francisco 2 ' 4 , P. J. S. Gil 2 ' 5 , and J. Paramos 1 - 6 
Instituto Superior Tecnico, 
Avenida Rovisco Pais 1, 1049-001 Lisboa, Portugal 
(Dated: October 20, 2008) 

We present a methodology based on point-like Lambertian sources that enables one to perform 
a reliable and comprehensive estimate of the overall thermally induced acceleration of the Pioneer 
10 and 11 spacecraft. We show, by developing a sensitivity analysis of the several parameters of 
the model, that one may achieve a valuable insight on the possible thermal origin of the so-called 
Pioneer anomaly. 



PACS numbers: 07.87.+V, 24.10.Pa, 44.40.+a 



I. INTRODUCTION 

A. General Background 

The existence of an anomalous acceleration on the Pi- 
oneer 10 and 11 spacecrafts, sun-bound and with a mag- 
nitude of api ~ (8.5 ±1.3) x 1CT 10 m/s 2 has been put 
forward a decade ago, using two independent code anal- 
yses [1, 2}. Attempts to account for these phenomena as 
a result of a misestimation of the systematic effects of 
thermal nature were first considered in Ref. [3, 4]. Pos- 
sible additional contributions, ranging from electric or 
magnetic forces, to mechanical effects or errors in the 
Doppler tracking algorithms used, have all be shown to 
be unsuccessful. 

Although initially dismissed, a much touted hypoth- 
esis for a physical explanation of the effect lies in the 
reaction force due to thermal radiation arising from the 
main bus compartment and the radiothermal generators 
(RTGs), either directly pointing away from the sun, or 
reflected by the main antenna dish. Clearly, an accelera- 
tion arising from the thermal dissipation should present a 
similar secular trend as the RTGs available power decay; 
regarding this point, one must note that another analysis 
has shown that such a signature in the anomaly may be 
found (i.e., it also possesses statistical significance), char- 
acterized by a linear decay with a time constant larger 
than 50 years [2]: given the ~ 88 years half- life of the 
plutonium source in the radio-thermal generators, which 
should be somewhat lowered due to degradation of the 
thermal coupling, this still leaves room for thermal radi- 
ation to account for the Pioneer anomaly. The latter is 
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being thoroughly examined by groups within the Pioneer 
collaboration team [5, 6]. 

In what concerns other effects, one can safely disregard 
electromagnetic forces, solar radiation and solar wind 
pressure as the cause for the anomalous acceleration [1]. 
Other sources for anomalous effects have been discarded, 
including the possibility that the Kuiper Belt's gravita- 
tional pull may give rise to the reported acceleration; this 
would require an abnormally high mass for this extended 
object, about two order of magnitude higher than the 
commonly accepted value of Micuiper = 0.3A/Earth [1, 7, 8] 
(for a variety of mass distribution models [7]). 

The two Pioneer probes are following approximately 
opposite hyperbolic trajectories away from the Solar Sys- 
tem. The fact that the same anomaly was found indicates 
a common origin to both spacecraft. This prompts for an 
intriguing question: what is the fundamental, and possi- 
bly new, physics behind this anomaly? 

Many proposals have been advanced to explain the 
anomaly as a previously undiscovered effect of new 
physics (see Ref. [9] and references therein, and also Refs. 
[10-12]). However, before one seriously considers the pos- 
sibility for new physics, an unambiguous description of 
the anomaly should be given. Unfortunately, the dis- 
tances at which the originally available Doppler measure- 
ments were conducted do not allow for a clear discrimina- 
tion of the direction of the acceleration: in particular, it 
is still not possible to discern between an acceleration to- 
wards the sun or the Earth, along the line of sight. Ascer- 
taining this would provide a relevant insight concerning 
the origin of the anomaly: a line of action pointing to- 
wards the sun would indicate a gravitational origin (since 
solar radiation pressure is manifestly too low to account 
for the effect), while a Earth-bound anomaly would hint 
at either a modified Doppler effect (due to new physics 
affecting light propagation and causing an effective blue 
shift) or an incorrect modeling of Doppler data, possibly 
due to mismodeled Earth orientation parameters, incor- 
rect ephemerides estimates, Deep Space Network (DSN) 
and software clock drifts, i.e., an unaccounted systematic 
effect. An intriguing possibility could be a "congenital" 
relationship between the Pioneer anomaly and the so- 
called flyby anomaly [13]. The anomalous acceleration 
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may also lie along the spin axis of the spacecraft: this 
would indicate that onboard, underestimated systematic 
effects to be held responsible for it; finally, an anomaly 
along the velocity vector would hint at some sort of drag 
effect. 

Regarding the latter, it is worth stating that this ad- 
ditional drag does not seem to be due to dark matter 
or dust distribution, since these are currently known to 
a good accuracy, and yield much lower effects. Con- 
versely, one may ask what density should the environ- 
ment have, so that a v 2 dependent drag force would 
account for the anomaly: a straightforward calculation 
shows that this should be of order 10~ 19 g/cm 3 (see, 
e.g., Ref. [7]); for comparison, the density of interplan- 
etary dust, arising from hot-wind plasma [14], is below 
10~ 24 g/cm 3 ; the density of interstellar dust (directly 
measured by the Ulysses spacecraft) is even smaller, at 
about 3 x 10~ 26 g/cm 3 . Also, a modification of geodeti- 
cal motion, hinting at an extension of General Relativity, 
could also account for a velocity dependent anomalous ac- 
celeration (see, e.g., Ref. [15] for a detailed discussion). 

Furthermore, it is clear that a careful study of secular 
and spatial trends should be carried out, aiming to re- 
late with possible thermal or engineering causes for the 
anomalous acceleration. The previously available data is 
likely to refer to an insufficiently long mission timespan, 
which does not allow for a clear discrimination of a hypo- 
thetical variation of the anomaly; to overcome this diffi- 
culty, recently recovered data of the full mission is being 
analyzed by distinct groups within the Pioneer collabo- 
ration team, with several approaches aiming to obtain 
convergent answers to the above questions (see e.g. Ref. 

Although initially disregarded, the issue of the Pio- 
neer anomaly has grown in and number of peer-reviewed 
publications, reflecting the increasing concern of the 
physics community. The characterization of any addi- 
tional anomalous acceleration was part of the scientific 
objectives of several mission proposals put forward to 
the recent ESA Cosmic Vision 2015-2025 programme 
[f 5, 17]; unfortunately, these efforts were ill-fated, leaving 
the community without the means to get a direct answer 
to this intriguing enigma. 



B. Previous Work 

A clear assessment of several systematic contributions 
to the overall acceleration may be found on Table I, ex- 
tracted from Ref. [1] . These baseline figures give a good 
measure of the different orders of magnitude of the vari- 
ous effects involved, and show that they do not account 
for the reported anomaly. As it turns out, unaccounted 
thermal effects are the most conspicuous sources of a sys- 
tematic effect. In Refs. [3, 4], estimates were performed 
for the heat dissipation of several spacecraft components, 
and claimed that a combination of several sources could 
account for the anomalous acceleration. In order to as- 



certain or disprove these and other claims, a more re- 
cent and thorough study has carried out the convoluted 
task of carefully modeling the Pioneer probes, in order 
to reproduce all relevant thermal effects with a sufficient 
accuracy [6]; a similar, independent effort is being under- 
taken by other groups within the Pioneer collaboration 
team. 

Although still preliminary, these attempts seem to in- 
dicate that thermal effects may account for up to one 
third of the total magnitude of the reported anomaly [18]. 
As we shall see this result is consistent with our own es- 
timates which indicate that thermal effects can account 
from about 35% to 67% of the anomalous acceleration. 
However, it is the authors' opinion that the many param- 
eter estimation and modeling strategies available up to 
now somehow cloud the overall picture, with the phys- 
ical significance being hindered by the technical depth 
of the thermal behavior reconstitution. For this reason, 
the present work attempts to drift somewhat away from 
the full modeling of every engineering detail, and directs 
its attention to the physical basis of the aforementioned 
thermal behavior. This stated, it is clear that our ap- 
proach is a complementary tool to the current endeav- 
ors: indeed, while a poorer modeling of specific details 
will reduce the overall confidence of the obtained results, 
the added simplicity, computational clarity and speed al- 
low for a convenient and much needed sensitivity analysis 
of the several relevant parameters. 

In this paper, we present the main features and the 
first results of a method based on point-like Lambertian 
sources. As we shall see, the presented method is already 
compatible with previous studies; further developments 
shall focus on a more detailed analysis of the reflectiv- 
ity effects, while still aiming at a good balance of model 
simplicity, computational speed and physical realism. 



II. SOURCE DISTRIBUTION METHOD 
A. Motivation and Rationale 

As discussed in the previous section, no definitive 
statements about the origin of the anomaly can be put 
forward until its full characterization. This justifies an 
intensive effort to recover and analyze the full flight data, 
and to develop approaches to understand the overall ther- 
mal behavior of the Pioneer probes, so to measure any 
previously unaccounted thermal radiation effects and to 
isolate, rule out, or constrain possibly remaining, yet un- 
known, effects. 

However, the authors feel that this pursuit should be 
countered with an approach focusing on the physical ef- 
fects directly relevant to the understanding of the prob- 
lem. The central issue is how thermal radiation is emit- 
ted, and reabsorbed and/or reflected, by the external sur- 
faces of the spacecraft and what is the resultant reaction 
force. Hence, instead of a complex finite elements model, 
that requires modeling of the whole spacecraft, we pro- 
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TABLE I: Error budget for the Pioneer 10 and 11, taken from Ref. [1]. 



Item Description of error budget constituents 



Bias Uncertainty 
10~ 8 cm/s 2 10~ 8 cm/s 2 



1 Systematics generated external to the spacecraft: 

a) Solar radiation pressure and mass 

b) Solar wind 

c) Solar corona 

d) Electro-magnetic Lorentz forces 

e) Influence of the Kuiper belt's gravity 

f) Influence of the Earth orientation 

g) Mechanical and phase stability of DSN antennae 

h) Phase stability and clocks 

i) DSN station location 

j) Troposphere and ionosphere 

2 On-board generated systematics: 

a) Radio beam reaction force 

b) RTG heat reflected off the craft 

c) Differential emissivity of the RTGs 

d) Non-isotropic radiative cooling of the spacecraft 

e) Expelled Helium produced within the RTGs 

f ) Gas leakage 

g) Variation between spacecraft determinations 

3 Computational systematics: 

a) Numerical stability of least-squares estimation 

b) Accuracy of consistency/model tests 

c) Mismodeling of maneuvers 

d) Mismodeling of the solar corona 

e) Annual/diurnal terms 



+0.03 



+ 1.10 

-0.55 



+0.15 



+0.17 



±0.01 
± < 10~ 5 
±0.02 
± < 10~ 4 
±0.03 
±0.001 
± < 0.001 
± < 0.001 
± < 10~ 5 
± < 0.001 



±0.11 
±0.55 
±0.85 
±0.48 
±0.16 
±0.56 
±0.17 



±0.02 
±0.13 
±0.01 
±0.02 
±0.32 



Estimate of total bias/error 



-0.90 



±1.33 



pose to develop a faster, more versatile approach based on 
a distribution of a few point-like thermal sources, simu- 
lating the thermal radiation emitted from the spacecraft, 
and analyzing the effect of radiation when emitted di- 
rectly to space or when reflected or absorbed by another 
surface of the spacecraft. This approach is complemen- 
tary to the ones based on finite element analyses and 
does not focus on the inner behavior of each component 
or surface, but instead attempt to isolate different con- 
tributions from the major constituents of the vehicles, 
namely the RTGs, antenna dish, and main bus compart- 
ment. 

There are several arguments justifying the interest and 
the effectiveness of the present approach. It is impossible 
to model the Pioneer spacecraft in a very precise way: 
it was built decades ago, accuracy of the blueprints or 



existing models is limited and the precise properties and 
degradation or damage of the materials, after decades in 
space, is unknown. This implies that even in the case of 
a full model of the spacecraft educated guesses will have 
to be done, limiting the accuracy of the obtained results. 

The impossibility of reliably describing several key pa- 
rameters should also limit the accuracy of any conclusions 
derived from a more encompassing approach. Specifi- 
cally, the limited temperature data (provided only by six 
sensors on the main bus and two sensors on the RTGs) 
and poor knowledge of the optical properties of the ma- 
terials introduce substantial uncertainties in the final re- 
sult, whatever is the adopted strategy. Thus, it is clear 
that the total and electrical power, which are well known, 
must be the fundamental parameters for any analysis. As 
will be shown, our approach is based on this principle. 
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Moreover, thermal radiation possibly contributing to 
the anomalous acceleration depends on the external sur- 
faces of the spacecraft and how the total power (and tem- 
perature profile) is distributed among them. The insula- 
tion of the spacecraft walls should limit the gradient of 
the temperature along each of the main external surfaces 
(except in special places as, e.g., the louvers, that can be 
modeled as separate sources if required) and make all the 
modeling of the details of the compartment unessential 
to address the problem. We argue that small details and 
small gradients in temperature of the spacecraft external 
surfaces will not affect the results considerably since, as 
we will see, the results are not too much affected by the 
number of point-like sources representing an extended 
surface (keeping the power constant). A small number 
of point-like sources can then be used to simulate any 
foreseen temperature gradient along each surface or a 
small localized extra source of radiation. From the ther- 
mal radiation point of view this is similar to a unevenly 
distribution of power by the few point-like sources rep- 
resenting the surface. Sensitivity analysis regarding the 
details of the spacecraft: shape modeling, temperature 
gradients, and total power emitted by each surface, can 
be then performed by varying the power assigned to each 
individual source in a prescribed way. 

The Pioneer spacecraft is spin stabilized and any reac- 
tion force component due to radiation will cancel away 
over time except in the direction of the axis of rotation. 
Most of the small contributions possibly not taken into 
account should be irrelevant since, due to the geometry 
of the spacecraft, most of them are expected to be nor- 
mal to the axis of rotation. This effect can be verified 
through the sensitivity analysis if slightly different radi- 
ation distributions by the sources lead to similar values 
of the anomalous acceleration, as expected. It should be 
noted that, as we are modeling relatively large radiating 
surfaces as point-like sources, the present model cannot 
provide too reliable information about the total reaction 
torque induced by the thermal radiation into the space- 
craft. 

It is clear that any study of this scope necessarily in- 
volves a large number of assumptions and hypotheses. 
Therefore, it is important to have the ability to quickly 
test a wide variety of scenarios and reach unambiguous 
conclusions about their plausibility: this sensitivity anal- 
ysis is crucially facilitated by the short computation time 
of the present method. In addition, the simplicity of the 
formulation keeps the involved physics visible throughout 
the entire process, allowing for scrutiny of every step. Fi- 
nally, we emphasize that the key goal of our effort is to 
perform a wide spectrum study of the parameter space 
for several physical properties relevant to the thermal 
modeling of the Pioneer probes. Our approach, while 
less comprehensive than a finite element model, allows 
for a direct interpretation of results, easy adaptability, 
as well as rather short computation times. 

Obviously, this endeavor would be incomplete if its self- 
consistency could not be assessed. Thus, before tack- 



ling the more interesting, physical case of the Pioneer 
anomaly, a set of test cases is performed to ascertain the 
effectiveness of the method. Furthermore, the choice for 
a point-like source approach should also be verified; this 
may be achieved by increasing the number of sources and 
observing the convergence of the relevant quantities and 
results. If deemed satisfactory, one may safely assume 
that continuous surfaces and components can be suit- 
ably modeled by point-like sources, so to still reproduce 
the physical interplay between them, and hence allow for 
an extrapolation to the Pioneer vehicles. 



B. Physical Formulation 

Our method is based on a distribution of isotropic and 
Lambertian point-like sources. If W is the emitted power, 
the time-averaged Poynting vector-field for an isotropic 
source located at (xo,yo, z o) is given by 



w 



(x-x ,y-yo,z- z ) 



4tt \( t . _ 
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(1) 



In the case of a Lambertian source the intensity of the 
radiation is proportional to the cosine of the angle with 
the normal 
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(2) 

Typically, one uses isotropic sources to model point-like 
emitters and Lambertian sources to model surfaces. The 
Poynting vector field of the source distribution is, then, 
integrated over the surfaces in order to obtain the amount 
of energy illuminating these, and the force produced. The 
former is given by the time-averaged Poynting vector flux 



£di um = / S • n dA = 



(3) 



/ S(G(M))-(- 



where the function G(s, t) parameterizes the relevant sur- 
face. The radiation illuminating the surface will produce 
a perpendicular force; integrating this force, i.e., the ra- 
diation pressure multiplied by the unitary normal vector, 
will give us the total force acting upon the corresponding 
surface. The radiation pressure is thus given by 



Prad = — S • n , 
c 



(4) 



taking into account a radiation pressure coefficient 1 < 
a < 2. The case a = 1 corresponds to full absorption 
while a — 2 indicates full diffusive reflection. 
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There will also be a force acting on the source of the 
radiation; this can be obtained by integrating the radi- 
ation pressure multiplied by a normalized radial vector 
field along a generic surface 

/"S-n S , „ 
F m * = J— m dA . (5) 

If the object in study has a reasonably complex geome- 
try (such as the Pioneer spacecraft) there will be shadows 
cast by the surfaces that absorb and reflect the radiation. 
The shadowing effect of the illuminated surfaces is cal- 
culated with this same expression and then subtracted 
to the force obtained for the emitting surface. Alterna- 
tively, one may use an integration surface that encom- 
passes the illuminated surfaces. The total result is the 
sum of all effects F; — force on the emitting surface, 
shadows and radiation pressure on the illuminated sur- 
faces a T h = J2i Fi/m Pio . 



C. Test Cases 

In order to demonstrate the efficiency of the proposed 
method, we define a set of test cases to assess the quality 
of our approximation. The key point is the ability to ad- 
equately represent the thermal radiation emitted from an 
extended surface by a small number of point-like sources, 
as opposed to having many small thermal radiating ele- 
ments. 

In the performed test cases, a square emitting surface 
with 1 m 2 is considered. The three components of the 
force are then computed: force on the emitting surface, 
shadow caused by another surface at a given position, and 
radiation pressure on the surface. We compare the results 
for different numbers of sources, while maintaining the 
total power fixed. It is expected that the result converges 
to the exact solution as the number of radiation sources 
increases. Our study shows that one is able to get a 
reliable error estimate even when using a small number 
of sources to model a surface. 

For a surface emitting radiation that does not illumi- 
nate other surfaces, one finds that the force is perpendic- 
ular to the former and only depends on the total emitted 
power. Using Eq. (5) with Lambertian sources on a sur- 
face on the Oxy plane, one obtains a force in the z-axis 
direction, of magnitude (2/3)W sur f/c. 

Computation of the shadow and pressure radiation on 
other surfaces yields results that are not independent 
from the source distribution. In order to acquire some 
sensibility on that dependence, we plot the variation of 
the radiation intensity with the elevation and the azimuth 
for 1, 4, 16, 64 and 144 source meshes, as depicted in Figs. 
1 and 2. 

A visual inspection of the results indicate that, even for 
one source, the maximum deviation occurs at the higher 
angles of elevation and is less than 10%. For the relevant 
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FIG. 1: Polar plot of the intensity variation with elevation 
of the radiation emitted by a surface on the Oxy plane (solid 
angle Q), when considering 1, 4 and 16 Lambertian sources 
(full, dashed and grey curves, respectively), maintaining the 
total power emitted by the surface constant at 1 IF (the curves 
for 64 or 144 sources overlap the one for 16 sources). The 
intensity at higher elevations (close to vertical) diminishes 
with the number of sources, compensating the slight increase 
at the lower angles. 



I(Wrad _1 ) 




FIG. 2: Same as Fig. 1 but for intensity variation with az- 
imuthal angle 6. All lines are superimposed, confirming that 
the total power is kept constant. 



angles for the Pioneer spacecraft configuration, devia- 
tions will be considerably smaller. In order to confirm 
this estimate, the force acting on a second 1 m 2 surface 
for several different positions is computed. A total of 
nine representative configurations were considered, with 
different positions and tilt angles, as summarized in Ta- 
ble II. The deviation between the 1, 4, 16, 64 and 144 
source meshes is then verified. 

TABLE II: Positions considered for the second surface in test 
cases. The first (emitting) surface is in the x—y plane centered 
at the origin. Considered distances between both surfaces are 
typical for the Pioneer spacecraft. 

Test case Surface center position Surface tilt angle 



# 


(m) 


(°) 


1 


(2, 0, 0.5) 


90 


2 


(2, 0, 1.5) 





3 


(2, 0, 1.5) 


30 


4 


(2, 0, 1.5) 


60 


5 


(2, 0, 1.5) 


90 


6 


(1, 0, 2) 





7 


(1, 0, 2) 


30 


8 


(1, 0, 2) 


60 


9 


(1, 0, 2) 


90 



Our study shows that the highest deviation occurs for 
Test Case 8, which confirms our expectation, since the 
second surface is set at high elevation from the emitting 
surface, as depicted in Fig. 3. The results in Table III 
show a difference of approximately 6% between the force 
obtained with one source and the results for the finer 
meshes (16, 64 and 144 sources). Nevertheless, the latter 
are all within 0.5% of each other, and the intermediate 4 
source mesh has a deviation of only 1.5%. 

TABLE III: Results for Test Case 8 (cf. Table II) considering a 
total emission of 1 kW. As the number of sources to represent 
the thermal emission of a surface change, the resultant force 
components appearing by shadow on the secondary surface 
remain almost the same. 



Sources 


Energy flux 


Force components (x, y, z) 


# 


(W) 


(10" 7 N) 


1 


45.53 


(2.016, 0, 2.083) 


4 


45.53 


(1.918, 0, 2.003) 


16 


45.53 


(1.895, 0, 1.984) 


64 


45.53 


(1.890, 0, 1.979) 


144 


45.53 


(1.889, 0, 1.978) 



For the typical angles of the Pioneer probe's configu- 
ration, one may take as figure of merit Test Cases 1 and 
3. For the first case, depicted in Fig. 4, the radiation 
pressure and shadow yield the results shown in Table IV. 
The analysis of these results shows that, for 16, 64 and 
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FIG. 3: Geometry of Test Case 8 (cf. Table II): thermal 
emission from a surface is simulated by a different number of 
Lambertian sources evenly distributed on the surface, main- 
taining the total power emitted constant, and the effect on 
a second surface is observed. This is the test case where the 
highest variation with the number os sources considered were 




FIG. 4: Same as Fig. 3, for Test Case 1. 



144 sources, the variation in the energy flux and force is 
less than 0.5%. In addition to that, the difference when 
compared with the results from finer meshes is less than 
5% for 1 source and less than 1.5% for a 4 source mesh. 
The results in Table V show, for Test Case 3, a variation 
of less than 5% between the results for 1 source and 144 
sources. The convergence is, as in both previous cases, 
achieved for the 16, 64 and 144 source meshes, with a 
variation of less than 0.25%. 

For all test cases examined, the convergence of the re- 
sults occurs at a similar pace and yields, for all cases, 
similar small deviations. Ultimately, we conclude that a 
4 source mesh, with deviations around 1.5%, would be 
adequate for the desired balance between precision and 
simplicity. These results provide a fairly good illustra- 



TABLE IV: Same as Table III, for Test Case 1. 



Sources 


Energy flux 


Force components (x, y, z) 


# 


(W) 


(10~ 7 N) 


1 


15.34 


(0.9300, 0, 0.1514) 


4 


15.92 


(1.028, 0, 0.1638) 


16 


16.09 


(1.038, 0, 0.1675) 


64 


16.13 


(1.040, 0, 0.1684) 


144 


16.14 


(1.041, 0, 0.1686) 


TABLE V: Same as 


Table III, for Test Case 2. 


Sources 


Energy flux 


Force components (x, y, z) 


# 


(W) 


/in -7 ivh 

In ) 


1 


19.20 


(0.4952, 0, 1.037) 


4 


19.83 


(0.5032, 0, 1.082) 


16 


19.99 


(0.5050, 0, 1.093) 


64 


20.03 


(0.5054, 0, 1.096) 


144 


20.04 


(0.5055, 0, 1.096) 




FIG. 5: Pioneer spacecraft model geometry considered in cal- 
culations, back view: high gain antenna and hexagonal main 
bus compartment. 



tion of the power of our method and how well we can 
estimate the radiation effects on the Pioneer probes. In 
particular, the deviation is always well below 10%, even 
with the roughest simplifications allowed by the chosen 
method. One may then conclude that, for the scales and 
geometry involved in the Pioneer anomaly problem, the 
source distribution method is, not only consistent and 
convergent, but that it provides a very satisfactory es- 
timate of the thermal radiation effects, even considering 
all uncertainties involved. 

Finally, after analysing the convergence of the method, 
we have also considered two additional test cases to assess 
the effect of ignoring some surface features, such as the 
equipment attached to the external walls of the space- 
craft. These results indicate that, unless large tempera- 
ture gradients are present, no significant errors will arise 
from considering flat surfaces and not taking into account 
all the details of the spacecraft. 



III. THERMAL RADIATION MODEL OF THE 
PIONEER SPACECRAFT 

A. Geometry 

The problem of modeling the Pioneer spacecraft can 
be considerably simplified with some sensible hypotheses. 
The first and most important one rests upon the fact that 
the probes are spin stabilized. Since it is also assumed 
that the probe is in a steady-state thermal equilibrium 
through out most of their journey, the time-averaged ra- 
dial components of any force generated by anisotropic ra- 
diation will be negligible. In addition, the probe's axis of 



rotation (taken as the z-axis) is approximately pointing 
towards Earth, which is also the approximate direction 
of the anomalous acceleration. 

In this study, we consider a simplified version of the 
spacecraft geometry, which retains only its most impor- 
tant features, as depicted in Figs. 5 and 6. Our model 
considers the RTGs, a prismatic equipment compartment 
and the antenna — a paraboloid, parametrized by the 
function G(s, t) = (s, t, a(s 2 + £ 2 )), with a parabolic 
coefficient a = 0.25 m -1 (c.f. Eq. (3)). Dimensions 
are taken from the available Pioneer technical drawings. 
Our results are obtained through the integration of the 
emissions of the RTG and lateral walls of the equipment 
compartment along the visible portion of the antenna. 
Note that radiation emitted from the front surface of the 
Pioneer cannot be reflected by other surfaces and will 
be counted as a whole. The surface of the compartment 
facing the antenna will be discarded for now as its con- 
tribution is fairly small for obvious geometric reasons: 
escaping radiation will be attenuated by multiple reflec- 
tions between these two components and will be mainly 
in the radial direction, not contributing significantly to 
the anomalous acceleration. The antenna itself is ex- 
pected to have a very low temperature (~ 70 K) with 
an approximately uniform distribution, not only axially, 
but also considering the front and back surfaces of the 
paraboloid (as visible in the results from Ref. [18]); there- 
fore, its contribution can be regarded as negligible, with 
the surface acting solely as a reflector for the incoming 
radiation. 

As we shall see, this simplified model captures the most 
important contributions to the thermal reaction force. 
The RTGs and the main equipment compartment are 
actually responsible for most of the emitted thermal ra- 




FIG. 6: Schematics of our simplified model of the Pioneer spacecraft, with relevant dimensions (in cm); second RTG and truss 
are not represented to scale, for convenience. Lateral view indicates the relative position of the RTGs, box compartment and 
the gap between the latter and the high gain antenna. 



diation. In the case of the equipment compartment, the 
most important contribution is from the louvers located 
in the front wall (facing away from the sun) with conse- 
quences for the total power distribution. 



B. Point-like Source Distribution 

In order to estimate the thermal effects, a separate 
analysis of the three main contributions must be per- 
formed. The front wall of the probe, where the louvers 
arc located, is perpendicular to the axis of rotation: its 
contribution corresponds to a force (2/3)Wf ron t/c point- 
ing in the sun-ward direction along the probe's axis. The 
contribution from the side walls of the main compart- 
ment is obtained from the integration of the shadow and 
radiation pressure components along the antenna. The 
shadow of the RTGs was neglected since they are small, 
relatively distant, and most of its effect would be in the 
radial direction. Following an approach similar to the 
one used in the test cases, in order to verify the conver- 
gence of the result, the integration is performed for an 
increasing number of sources. The results converge fairly 
quickly and the deviations are all below 2.5%, confirm- 
ing the consistency previously demonstrated in our test 
cases. The obtained values show that between 16.8% to 
17.3% of the power emitted from the side walls of the 
compartment is converted into a sun-ward thrust along 
the z-axis. 

It is also important to verify how the results are af- 
fected by a non-uniform temperature distribution. This 
is simulated by varying the relative power of the point- 
like sources in each surface, keeping the total power at- 



tributed to the surface constant. A variation of 20% 
in power between sources (simulating a 5% temperature 
variation) gives no significant changes in the final result 
— with relative differences smaller than 1%. 

Finally, the RTG contribution is computed through 
two different models. The first, simpler scenario, mimics 
each RTG with a single isotropic source. In this case, 
the point- like source has the whole power of the RTG. 
In the second model, the cylindric shape of the RTG is 
taken into account and a Lambertian source is placed at 
each base. Actually, it is only necessary to consider the 
source facing towards the centre of the spacecraft, as the 
remaining RTG radiation will be emitted radially and 
its time-averaged contribution vanishes. In this case, the 
Lambertian source has a certain amount of the total RTG 
power, as discussed in the following sections. Depending 
on the model considered, either 1.9% of the total power or 
12.7% of the power emitted from the base of the cylinder 
(equivalent to approx. 2% of total RTG power, if the 
temperature is uniform) is converted into thrust. 

These preliminary results do not take into account dif- 
fusive reflection, as allowed by Eq. 4. In the subsequent 
section, more accurate results will be presented and dis- 
cussed. 



C. Available Power 

The available power on the Pioneer spacecraft is one 
of the few measured or inferred parameters that is rea- 
sonably well known. In addition, it is physically more 
consistent to consider the power instead of the tempera- 
ture readings as it is the independent variable from which 
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all estimates of the resulting thermal effects are derived. 
Of course, the energy balance to the spacecraft in steady- 
state conditions relates the temperature Tj of a surface i 
with the power budget of the spacecraft 



-^absorb + E gcn — AitiaT^ 



(6) 



where Ai are the relevant areas and e, the emissivity of 
each surface i. 

Notice that, since the optical properties of the surfaces, 
as well as their evolution with time, are not well known, 
temperature estimates are quite uncertain. Furthermore, 
as a variation in the emissivity implies the violation of 
the conservation of energy, a new solution for the tem- 
perature must be obtained iteratively for each different 
set of optical properties used so to maintain the correct 
power. 

All the power generated onboard the probes comes 
from the two plutonium-238 RTGs. Given that just a 
fraction of the generated heat is converted into electric- 
ity, the remaining power is dissipated as heat. There 
will be some conduction of this heat through the truss to 
the central compartment, however, considering the small 
section of this structure, it is reasonable to admit that 
it will have a reduced impact on the total RTG radiated 
power [19]. It is, thus, considered that all of the RTG 
thermal power is dissipated as radiation from the RTG 
itself. 

The electrical power is consumed by the various in- 
struments located in the main compartment, despite a 
considerable portion of it being used in radio transmis- 
sions from the high gain antenna. As mentioned in Ref. 
[1] , the total RTG thermal power at launch was 2580 W, 
producing 160 W of electrical power. This means that, 
at launch, approximately 2420 W of thermal power has 
been dissipated by the RTGs. Taking into account the 
plutonium decay with a half-life of 87.74 years, the total 
on-board power variation with time (in years) is given by 



TF tot (t) = 2580exp(-J^|) 



(7) 



IV. RESULTS AND DISCUSSION 
A. Order of Magnitude Analysis 

Before undertaking a more rigorous numerical esti- 
mate, one may use the results described above to per- 
form a preliminary order of magnitude analysis. This 
allows one to obtain a concrete figure of merit for the 
overall acceleration arising from thermal effects, which 
can be compared with the api ~ 10~ 9 m/s 2 scale of the 
Pioneer anomaly. 

From the spacecraft specifications, one has a total mass 
mpio ~ 230 kg, and separate RTG and equipment com- 
partment powers IFrtg ~ 2 kW and W e quip ~ 100 W. 
As already discussed, the integration of the emissions 
of the RTG and instrument compartment indicate the 
proportion of emitted power that is effectively converted 
into thrust. If we consider the simpler model discussed 
is section III B and the power emitted from each surface 
proportional to its area (equivalent to assuming uniform 
temperature and emissivity in the RTGs and equipment 
compartment), we obtain 



^rtg - 2 x 10 

^sides rs "' 10 
F front - 2 X 10 



> Wrtg 



(8) 



iWe, 



C 

. 1 W'equip 



One can easily estimate the acceleration of the space- 
craft due to the thermal effects arising from the power 
dissipation of the RTGs and equipment compartment: 



., , n _2 " :n <; ,. m - n> 
ortg ~ 2 x 10 ... i.. 



3 x 10" 



6 x 10" 1U m/s 2 , (9) 

mpi c 

1 Wequip m-10 



m Pio c 



4.4 x 10~ iu m/s 2 . 



This clearly indicates that both contributions are rele- 
vant to account for the reported anomalous acceleration 
of the Pioneer probes. Furthermore, it also shows that 
the RTGs and the instrument compartment yield similar 
thermal effects, so that one cannot focus solely on one 
of these sources when modeling the spacecraft (this had 
already been revealed by the analysis of Ref. [6]). 



Telemetry data reveals that the electrical power decays 
at a faster rate than the plutonium radioactive decay; in 
the latest stages of the mission, about 65 W were avail- 
able. Most of the electrical power is dissipated inside 
the main compartment. The electrical heat in the space- 
craft body was around 120 W at launch, dropping to less 
than 60 W at the latest stages of the mission [6] , follow- 
ing an approximate exponential decay with a half-life of 
about 24 years. This difference in decay rates is mostly 
attributable to thermocouple degradation. 



B. Thermal Force Estimate 

Encouraged by the estimate outlined above, one may 
now proceed with a more thorough evaluation of the ex- 
isting thermal effects, using our point-like source model- 
ing. 

In this section we shall use a model with 4 sources in 
each side panel of the equipment compartment and Lam- 
bertian sources at the bases of the RTGs, as discussed in 
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section IIIB. We believe this model gives us the best 
compromise between accuracy and computation time — 
the deviation of the source distribution relative to the 
finer meshes is less than 0.5 %. Integrating the radiation 
pressure and shadow components using the methodology 
presented in section II B and extracting the axial com- 
ponent, we obtain an expression that yields the thermal 
acceleration 

(0.168W/ sidcs + |W front + 0.128T^ RTG b) , in , 

a Th = ^ 2 L j (io) 

m Pio c 

where W^ides and Wf r0 nt are the powers emitted from the 
side panels and front of the equipment compartment and 
WftTGb is the power emitted from the base of the RTG 
facing the centre of the spacecraft. Remaining contribu- 
tions are much smaller, as discussed in sections III A and 
IIIB. 

A critical analysis of this expression, bearing in mind 
the spacecraft geometry, reveals that all considered con- 
tributions yield a sun-ward acceleration: the Upfront com- 
ponent radiates directly in a direction away from the sun, 
while the other two components W s idcs and W^RTGb radi- 
ate laterally, illuminating the high gain antenna — which 
will yield a significant shadow and radiation pressure. 
The question now resides in correctly estimating each one 
of these powers. We shall consider the 1998 readings, as 
found in the graph of Ref. [6] , namely: Wrtg = 2050 W 
and Woquip = 58 W. These are the dissipated thermal 
powers at the RTG and equipment compartment. 

The simplest scenario, with uniform temperature and 
optical properties (emissions proportional to the surface 
area, as in the previous section), leads to 



aih/10" 



Wsides - 21.75 W 
W front = 18.12 W 

WRTGb = 41.11 W 



(11) 



yielding an acceleration axh = 3.05 x 10~ 10 m/s 2 . This 
amounts to about 35% of the anomalous acceleration. 
However, it is wise to undertake a critical analysis of 
this figure: considering the available temperature maps 
of Rcfs. [6, 18], we see that the temperature anisotropics 
along the sides of the equipment compartment fall within 
the tested cases, as discussed in section IIIB. However, 
the RTG temperature distribution should deserve further 
attention, as there are significant temperature changes 
between the wall of the cylinder, the bases and the fins. 
In addition, it is expected that the front wall of the equip- 
ment compartment will have a larger contribution than 
the side walls, due to the presence of the louvers. 

Taking these considerations into account, one can ana- 
lyze the variation of the emitted power in the louvers and 
at the base of the RTG, since these are the two critical 
parameters in the calculation. If we consider that the 
louvers are closed and have a similar cmissivity to the 
rest of the equipment platform, we can plot the variation 




3.0 Tequip 



FIG. 7: Variation of the resulting acceleration with the tem- 
perature ratio between the louvers and the equipment plat- 
form, considering similar emissivities for both multi-layer in- 
sulations. 
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FIG. 8: Variation of the resulting acceleration with the tem- 
perature ratio between the base of the RTG cylinder and the 
fin temperature. 



of the acceleration with the temperature ratio between 
the louvers and the mean temperature of the platform, 
while keeping the total power constant. This is depicted 
in Fig. 7. One can perform a similar analysis for the 
RTGs, considering the ratio between the temperatures 
at the base of the cylinder and the fins (Fig. 8). 

Figs. 7 and 8 are illustrative of the main strength of our 
method: it allows for a fairly quick and accurate analysis 
of the dependence of the final result on different param- 
eters. Through Eq. (10) and sensible variation of the 
power parameters, one can match temperature readings 
and consider hypotheses for the variation of the optical 
properties. 

We can now perform a second estimate considering the 
RTG cylinder bases and wall as having a 15% and 30% 
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higher temperature than the fins, respectively. Assum- 
ing also that the closed louvers have similar emissivities, 
although a 100% higher temperature than the rest of the 
equipment compartment could be possible, one obtains 
the following values for the powers: 

Wsides = 9.97 W , (12) 
VKfront = 39.71 W , 
VK RTG b = 49.67 W . 

In this case, one can account for 57% of the anomalous 
acceleration, that is, axh = 5.00 x 10~ 10 m/s 2 . 

So far, our results do not consider reflections, i.e., full 
absorption of the radiation by the illuminated surfaces. 
In this study, we shall introduce diffusive reflection by 
assigning a value to the a parameter in Eq. (4). For the 
kind of aluminum used in the construction of the antenna, 
the reflectivity is, typically, around 80% for the relevant 
wavelengths, yielding a = 1.8. For the multi-layer insula- 
tion of the equipment platform, a value of a = 1.7 is con- 
sidered. In these conditions, the illumination factors in 
Eq. (10) are modified to account for the reflection. With 
the same temperature conditions as in the previous case, 
the resulting acceleration is axh = 5.75 x 10~ 10 m/s 2 - 
approximately two thirds of the anomalous acceleration. 

The results presented in this section give us a fairly 
good idea of the changes involved when considering dif- 
ferent hypotheses and parameters. The three discussed 
scenarios here illustrate how one can use our method to 
identify the most sensitive parameters and quickly assess 
the effect of the existing uncertainties, suggesting where 
models must be refined in order to increase confidence in 
results. 



V. CONCLUSIONS 

In this work we have developed a method to account for 
the acceleration of the Pioneer spacecraft due to thermal 
effects, based on point-like Lambertian sources. The flex- 
ibility and computation simplicity of our method allow 
for a reliable and fast estimate of the acceleration due to 
the various thermal contributions of the spacecraft com- 
ponents. This is sharply contrasting with the complexity 
and computationally demanding nature of the finite el- 
ement analysis. Our methodology is potentially useful 



for a thorough parametric study of the various thermal 
contributions, as discussed in sections III and IV. 

Our method allows for a reasonable degree of accu- 
racy and the numerical error estimates provided by the 
numerical calculation package are of the order of 10~ 14 
or less, while the approximation of the geometry with 
point-like sources results in a deviation of less than 1%, 
as discussed in sections II C and III B. This should not be 
understood as an indication of the accuracy of the result- 
ing accelerations, when compared to the reported case of 
the Pioneer anomaly, but as a measure of self-consistency 
and reliability of the developed method which should 
be expanded to model the physical system of the Pioneer 
spacecrafts more closely, while maintaining the desired 
flexibility and computational speed. 

We find, after identifying the main contributions for 
the power of the various components of the spacecraft 
(RTGs, antenna and equipment bus compartment), fig- 
ures ranging from 35% to 57% of the anomalous acceler- 
ation disregarding reflection. Inclusion of reflection im- 
plies that one can account for about 67% of the anomaly. 

The natural continuation of this work will involve the 
refinement of the geometrical modeling, including the 
specular component reflection. In addition, and possibly 
more relevantly, we aim to pursue the identification of 
parameters that most significantly affect the final result 
— namely temperatures, emissivities and reflectivities of 
the various components, such as the louvers and the RTG 
case. In any case, our analysis does achieve a reasonable 
level of agreement with other thermal models based on 
finite element methods [6, 18]. 
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